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Abstract. Jets can be used to probe the physical properties of the high energy density matter created in 
collisions at the Relativistic Heavy Ion Collider (RHIC). Measurements of strong suppression of inclusive 
hadron distributions and di-hadron correlations at high pr have already provided evidence for partonic 
energy loss. However, these measurements suffer from well-known geometric biases due to the competition 
of energy loss and fragmentation. These biases can be avoided if the jets are reconstructed independently 
of their fragmentation details - quenched or unquenched. In this paper, we discuss modern jet recon- 
struction algorithms (cone and sequential recombination) and their corresponding background subtraction 
techniques required by the high multiplicities of heavy ion collisions. We review recent results from the 
STAR experiment at RHIC on direct jet reconstruction in central Au+Au collisions at ^^snn = 200 GeV. 

PACS. 21.65.Qr Relativistic Heavy Ion Collisions - 13.87.Ce Jet Production 



1 Introduction 

A highly collimated "spray" of particles also known as 
"jets" is produced in high energy collisions. Hadronic jets 
' are the experimental signatures of quarks and gluons (par- 
' tons) and they reflect the underlying parton kinematics [1^ 
, [2l[3l[l]. Cross section measurements of jets are performed 
' at many hadronic and leptonic colliders to check in de- 
[ tail perturbative QCD (pQCD) calculations, to help de- 
termine parton distribution functions and to look for new 
physics. The inclusive jet cross section at Tevatron is mea- 
sured very precisely over 20 orders of magnitude and is 
found to be in a very good agreement with the NLO pQCD 
calculations using CTEQ 6.1 parton distribution functions 
[5l|6]. The robustness of the theoretical calculations on jet 
cross sections in p -\- p collisions motivates the use of jets 
as direct probes of partonic energy loss in dense matter 
generated in ultra-relativistic heavy ion collisions at RHIC 
and in near future at the LHC |7ii8j. 

Due to the limited pseudo-rapidity {r]) and azimuthal 
(j) coverage of the electromagnetic calorimeters of RHIC 
experiments during the first 6 years of RHIC operation, 
the background determination needed for direct jet re- 
construction was not possible. Instead other observables 
such as the strong suppression of inclusive hadron distri- 
butions and di-hadron correlations at high pt were mea- 
sured. However, such measurements suffer from geometric 
biases due to the competition of energy loss and fragmen- 
tation - the leading particle spectrum is dominated by 
relatively low energy jets that happen to lose little energy 
in the medium and fragment into higher pt particles [9j . 
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These biases can be removed if the reconstructed partonic 
kinematics is independent of whether the fragmentation 
is modified by the medium or not. An unbiased jet re- 
construction measurement in heavy ion collisions would 
give access to the full spectrum of fragmentation topolo- 
gies without geometric biases, enabling full exploration 
of quenching dynamics. In addition, fully reconstructed 
jets allow the measurement of qualitatively new observ- 
ables such as jet shapes, fragmentation functions, and en- 
ergy flow. Since 2006, the STAR barrel electromagnetic 
calorimeter (BEMC) with full azimuthal coverage and unit 
pseudorapidity acceptance is operational. This enables the 
study of the underlying event background required for full 
jet reconstruction in heavy ion collisions for the first time 
at RHIC. The experimental results discussed in this ar- 
ticle were presented for the first time during the Hard 
Probes 2008 meeting [lOj for the first direct measurement 
of jets and [llj for the accompanying measurement of jet 
fragmentation studies in heavy ion collisions. 



2 Jet Reconstruction Algorithms 

During the last 20 years, various jet reconstruction al- 
gorithms were developed to combine measured particles 
into jets. For a detailed overview of jet algorithms in high 
energy collisions, see p!2 l [T3 l[T4l [T5] and references therein. 
The primary requirement for jet algorithms is low sensitiv- 
ity to hadronization, radiation and splitting. They should 
also be defined equally at hadron and parton level. This 
was suggested by Weinberg 30 years ago [4 . According to 
him "Quark and gluon jets can be compared to detector 
jets, if jet algorithms respect collinear and infrared safety." 
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Full jet reconstruction in heavy ion collisions is a new 
frontier. As in leptonic and hadronic collisions, in heavy 
ion collisions the chosen jet reconstruction method should 
also be theoretically and experimentally consistent. As the 
experimental energy determination degrades the resolu- 
tion somewhat with the detectors, the algorithm should 
also aim to minimize resolution effects from unrelated 
sources such as the underlying event. The expected in- 
crease in Large Hadron Collider (LHC) luminosities (20 
to 200 collisions in a detector) requires that the tradi- 
tional jet algorithms have to be improved with underlying 
event subtraction techniques for p + p collisions to resolve 
events for pile up. These improved techniques can also 
be used in heavy ion environments (Au+Au or Pb+Pb) 
where the background subtraction is required due to large 
multiplicities of produced particles [14^16j. 

The algorithm should also be detector independent 
i.e., allow the combination of particles detected in var- 
ious detectors. The results discussed in this article are 
from jets that are reconstructed by combining the neu- 
tral energy from the BEMC and charged particles from 
the Time Projection Chamber (TPC) of the STAR exper- 
iment. STAR'S TPC and BEMC detectors cover full az- 
imuth (0 < < 27r) and mid rapidity (— 1 < < 1) of the 
events. The Figure [T] shows an example of a reconstructed 
di-jet event in the STAR experiment. 



2.1 Cone Algorithms 

The cone algorithms have been used as a primary tool 
to identify jets at hadron colliders since the early 1980s. 
This algorithm is based on the picture that a jet consists 
of a large amount of hadronic energy in a small angular 
region. Therefore, the main method is to combine particles 
in r] — (j) space with their neighbors within a cone of radius 
R{R= ./A^- 
in Figure [21 



Arp'). This is illustrated schematically 




Tl 



Fig. 2. The transverse momentum of tracks and/or deposition 
of energy in calorimeter towers are represented as the solid 
squares in r\ and space. Red squares defined as seeds are the 
ones above a given threshold. Jet cones are solid circles around 
the seeds. All the energy of the particles is added for the given 
cone around the seed particles to estimate the jet energy. The 
background in the jet energy is estimated with the average 
of the total energy in cones without seeds (shown in dashed 
circles) and is subtracted on an event- by- event basis. 



To optimize the search and effectiveness of jet finding, 
these algorithms use splitting, merging, and iteration steps 
in the events of leptonic and hadronic collisions. However, 
to avoid instabilities in cone-finding due to large heavy 
ion background, we use a simple seeded cone without it- 
eration or split-merging steps, with cone radius R = 0.4 
and minimum seed of 4.6 GeV. The choice of the rela- 
tively small cone size is to suppress the underlying heavy 
ion background [17^18j. In p -\- p collisions ~ 80% of the 
jet energy is observed to be within R~ 0.3 for 50 GeV jets 
in the Tevatron data [6j. However, broadening of the jet 
fragmentation due to quenching in the medium formed in 
heavy ion collisions may reduce the fraction of the mea- 
sured energy in a given cone size and needs further explo- 
ration. To reduce the heavy ion background, the minimum 
accepted transverse momentum of charged particles, and 
the transverse energy of the calorimeter cells (p^*) is var- 
Corrections for double-counting of energy due to hadronicied between 0.1 to 2 GeV. This threshold cut does not 
energy deposition in the BEMC and to electrons are ap- remove all the background contamination on the jet en- 
plied. Two kinds of jet reconstruction algorithms are util- ergy and additional subtraction is needed. As presented 
lized; seeded cone (leading order high seed cone (LOHSC)) in Figure [2] schematically, the residual background is cor- 
and sequential recombination (kx and Cambridge /Aachen), rected based on the out-of-cone energy for the same p^^^ 
In the following, we briefly discuss these two algorithms averaged over the STAR acceptance but measured on an 
and the corresponding underlying event subtractions. event-by-event basis, and scaled to the cone area. 




Fig. 1. 21 GeV di-jet reconstructed from a single event with 
a combined transverse momentum and energy for charged and 
neutral particles per grid cell in the 77 and plane from 0-20% 
most central Au+Au collisions 111 . 
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A recently developed seedless infrared-safe cone algo- 
rithm (SISCone) [19j resolves most of the ambiguities as- 
sociated with the various cone algorithms. The SISCone 
algorithm is already used in p -\- p collisions at = 200 
GeV and the first results can be found in 20 . 



2.2 Sequential Recombination Algorithms 

The sequential recombination algorithms have been used 
extensively in the Tevatron as they are collinear and in- 
frared safe [6l[2TJ[22j. In these types of algorithms, arbi- 
trarily shaped jets are allowed to follow the energy flow 
resulting in less bias on the reconstructed jet shape than 
with cone algorithms [14]. Figure [3] represents a schematic 
comparison of the jet areas for cone and kx type algo- 
rithms. 




kyjet Cone jet 



Fig. 3. A schematic comparison of the jet area for cone and 
kx type algorithms. 



The sequential recombination algorithms combine ob- 
jects in relative to the closeness of their px- Particles are 
merged into a new cluster via successive pair- wise recombi- 
nation. Algorithmic details can be found in \6^2T] and ref- 
erences therein. The Fast Jet code package for sequential 
recombination algorithms was used for the STAR analy- 
ses in p + p and Au + Au collisions ^lQ l[TT|[2Qp3^ [24] . This 
package includes kx, Cambridge/ Aachen (CAMB), anti- 
kx, and an interface to external jet finders such as Seed- 
less Infrared Safe Cone (SisCone) via a plugin mechanism 
[19]. For infrared and collinear safe algorithms an active 
area (Aj) of each jet is estimated by filling an event with 
many very soft particles and then counting how many are 
clustered into a given jet. If the underlying event is dis- 
tributed uniformly in rj and (j) then this noise density can 
be subtracted from the measured jet energy on an event- 
by-event basis to correct for the background energy under- 
lying the jet. In simulations, this correction is observed to 
recover most of the jet energy when they are reconstructed 
in pile up and heavy ion backgrounds [14j. The kx, Cam- 
bridge/Aachen and anti-kx algorithms are all based on the 
same sequential recombination algorithm, but they differ 
in the distance measure that is used to group particles in 
to jets. 

Figure [4] shows the jet-area from kx algorithm from 
Fast Jet code for 0-10% central Au+Au events (MB-Trig), 
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Fig. 4. Jet area from kx algorithm is reconstructed uti- 
lizing the Fast Jet code package for real jets in mini- 
mum bias triggered 0-10% central Au + Au collision (MB- 
trig), in PYTHIA isolated jet events embedded in real central 
Au-\-Au events (PyEmbed) and in PYTHIA isolated jet events 
(PyTrue). 



next section for details) [10]. PYTHIA jets embedded in 
real Au+Au background events are observed to have the 
same area as jets from real Au+Au events. See the solid 
blue and dashed red histograms in Figured! The reduction 
in the MB-Trig jet area relative to PyTrue is well under- 
stood for sequential recombination algorithms on theoret- 
ical grounds [M]- 



3 Jet Reconstruction Analysis 
3.1 Event Selection and Terminology 



This analysis utilizes events of Au+Au collisions at ^/snn - 
200 GeV recorded by STAR. Using the multiplicity mea- 
surement performed with the STAR TPC, only the most 
central (0-10%) Au+Au collisions are selected. Two event 
sets were analyzed, based on fast on-line trigger configu- 
rations: 

(i) MB-Trig, (minimum-bias trigger) utilizing the coin- 
cidence between the two calorimeters at beam rapidity 
(Zero Degree Calorimeters) with signals in each greater 
than ~ 40% of the most probable amplitude for a single 
neutron, and 

(ii) HT-Trig (high-tower trigger) that satisfies the MB- 
Trig conditions and the additional requirement of 2-tower 
EMC clusters having at least a 7.5 GeV energy deposition. 

Three million MB-Trig events, corresponding to 300 
thousand 0-10% most central events are used for this study. 
A total of 80 million MB-Trig events were recorded by 
STAR during the 2007 Au + Au run, but only the event 
set used have been fully reconstructed off-line. The HT- 



compared to generated PYTHIA jets (PyTrue) and PYTHIA Trig is designed to enhance the recorded rate of high pt 
jets embedded in heavy ion background (PyEmbed, see photons and electrons. It may also serve to enhance the 
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recorded rate of jets. The HT-Trig data-set has been re- 
constructed in its entirety, corresponding to 500 /ib~^. 

In order to assess jet reconstruction energy resolu- 
tion, background subtraction, efficiency and acceptance, 
Monte-Carlo model studies based on PYTHIA 8.107 [25] 
are performed. PYTHIA events with high Et jets are gen- 
erated in three different configurations: 

(i) PyTrue: PYTHIA isolated jets including all particles 
except neutrinos. Jets are reconstructed using the PYTHIA 
internal jet finder, PyCell, for the cone algorithm, and 
Fast Jet for the sequential recombination algorithm. 

(ii) PyDet: PYTHIA isolated jets (parameterized detector 
response level) reconstructed using the jet algorithms that 
are also applied to the real data. 

(iii) PyEmbed: PyDet that are embedded in a background 
of real Au+Au 0-10% central events, with jets reconstructed 
with the jet algorithms that are also applied to the real 
data. 

In all the real events and simulations, only the highest 
energy jet per event is selected as the reconstructed jet. 

3.2 Energy Resolution 



CO 




Fig. 5. Distributions showing energy resolution; black AE = 
EpyDet - EpyTrue, red AE = EpyEmbed " EpyTrue and greeu 

AE — EpyEmbed — EpyDet • 

though smaller in magnitude, are also observed for the kx 
and Cambridge /Aachen algorithms. 



The energy resolution for jet reconstruction with vari- 
ous algorithms has been studied with isolated jets simu- 
lated with PYTHIA [25[l26l l27j. Figure [5] shows the event 
by event comparison of PyTrue, PyDet and PyEmbed 
from LOHSC algorithm. See [10] for comparison of the 
energy resolution with the kx and Cambridge /Aachen al- 
gorithms. Applied cuts and jet energy are specified in the 
figures. A shift of median due to un- measured particles, 
primarily neutrons and and the applied cut (hence 
loss of jet energy), is observed for the AE = EpyDet — 
EpyTrue histogram. Background effects are simulated using 
PYTHIA jets that are embedded in real Au + Au events. 
The distribution in Figure [5] is convoluted with the true 
jet spectrum to produce the observed jet spectrum. The 
effect of the heavy ion background on the jet energy can 
be seen in the AE = EpyEmbed — EpyXrue distribution. A 
positive AE in this distribution can distort the measured 
inclusive jet spectrum substantially, increasing the appar- 
ent yield at high Et and resulting in a harder spectrum. 
A correction to the spectrum must be applied to account 
for this effect. 

The influence of energy resolution on jet spectrum can 
be observed with PyDet, PyEmbed and PyTrue distri- 
butions shown in Figure [6] for the LOHSC algorithm for 
pcwt = 0.1 GeV. A large difference between PyEmbed and 
PyDet is observed as expected from the tail at positive 
AE due to large background in the red distribution in 
Figure [51 When the px threshold is increased, the back- 
ground fluctuations are reduced and the enhancement in 
the spectrum relative to the case without background is 
reduced to a negligible level [10]. Jet reconstruction in 
0-10% most central Au+Au collisions is similar to that 
of p+p collisions with a larger GeV threshold re- 
quirement. However, a reduction in the measured jet en- 
ergy (p^^* dependent bias) is introduced. Similar effects. 




ET[GeV] 



Fig. 6. Inclusive jet spectrum for PyDet, PyEmbed and 
PyTrue using the LOHSC algorithm. The p^^' = 0.1 GeV 
on track momentum and calorimeter cell energy is applied for 
PyDet and PyEmbed. Note the lower threshold on generated 
jet energy E^^'^^^^ > 5 GeV, which affects the reconstructed 
spectrum up to Et = 20 GeV. 



3.3 Jet Spectra Corrections and Comparisons 

The correction factors for the jet spectrum are estimated 
using PYTHIA simulated jets embedded in real Au+Au 
collisions . The Et dependent ratio of PyEmbed to PyTrue 
is calculated from Figure [6l A polynomial function fit to 
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Table 1. Correction factors for the inclusive jet spectrum, 
for different reconstruction algorithms and values of Pj^*. The 
range of values given indicates the correction factor variation 
from the lowest to highest jet Et shown in the figures. 
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the ratio distribution is used as a multiplicative correction 
to the inclusive spectrum. Table [1] shows the inclusive jet 
spectrum correction factors for various values. For the 
sequential clustering algorithms, the correction factors are 
closer to unity. 

The corrected inclusive jet spectrum for the LOHSC 
algorithm for the p^^ = 1 GeV is presented in Figure [71 
The p^^ = 1 GeV is selected as it corresponds to cor- 
rection factors close to unity. At this p^^ the competing 
effects of energy loss due to momentum threshold cut and 
the kick up in the jet spectrum due to the positive tail of 
energy resolution cancel each other. The solid triangles are 
for the MB-Trig data set and are corrected for resolution, 
acceptance and efficiency. 
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Fig. 7. Jet yield per event vs Et for 0-10% central Au + Au 
collisions, compared to the distribution from p -\- p collisions 
scaled by Neinary [29]. Triangle symbols are from MB-Trig and 
corrected for efficiency, acceptance and energy resolution. Open 
circles are from HT-Trig and not corrected for trigger bias. 
Only statistical error bars are shown for the Au + Au data. 
SoUd black squares are the distribution from p -\- p collisions, 
scaled by N Binary The yellow band represents the systematic 
uncertainty of the p + p measurement. 



In order to assess the biases in the jet spectrum re- 
constructed in central Au+Au collisions, we compare to 
the spectrum measured in p+p collisions. To account for 



nuclear geometric effects we scale the p+p spectrum by 
Neinary, the uumbcr of binary nucleon+nucleon collisions 
equivalent to a central Au+Au collisions, as calculated by 
a Glauber model [28] . The cross section for hard processes 
is expected to scale with Neinary if no nuclear effects are 
present. In the case of jet reconstruction, Nsinary scaling 
is expected if the reconstruction is unbiased, i.e. the jet 
energy is recovered independent of the particular mode of 
fragmentation, even in the presence of strong jet quench- 
ing. The N Binary scalcd jet spectrum from p + p collisions 
is shown in solid squares [29j. The yellow band represents 
the systematic uncertainty of the p + p jet measurement. 
Heavy ion jet spectrum is observed to agree with N Binary 
scaled p+p measurement within the ~ 50% systematic 
uncertainty of the normalization. 

The open circles in Figure [71 show the uncorrected jet 
spectrum from HT-Trig data which is substantially lower 
than the corrected MB-Trig spectrum. The correction fac- 
tors of energy resolution, efficiency and acceptance for the 
HT-Trig are expected to be small. A large trigger bias due 
to the additional 7.5 GeV energy deposition in EMCAL 
requirement for HT-Trig relative to MB-Trig is seen to 
persist at least to 30 GeV. Further statistics of MB-Trig 
data is needed to assess the bias at high p^ . 
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Fig. 8. Jet yield per event vs Et for 0-10% central Au + Au 
collisions obtained by the kx algorithm. The distribution from 
p + p collisions are scaled by Nsinary [22] • Triangle symbols 
are from MB-Trig and corrected for efficiency, acceptance and 
energy resolution. Only statistical error bars are shown for the 
Au + Au data. Solid black squares are the distribution from 
p + p collisions, scaled by N Binary The systematic uncertainty 
of the p + p jet spectra normalization is ^ 50%. 



Figures I8|9|10l show the comparison of inclusive jet 
spectra from the MB-Trig Au + Au data and the N Binary 
scaled p+p for the pr = 0.1 GeV threshold cut for kT, 
Cambrdige/ Aachen and LOHSC algorithms. While the 
agreement between N Binary scaled p+p and MB-Trig mea- 
surement is good for p9^^ = 0.1 GeV, it is also seen to be 
poorer with the larger pt threshold cut [10] . This suggests 
that p9^^ introduces biases which are not fully corrected 
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Fig. 9. Jet yield per event vs Et for 0-10% central Au + Au 
collisions obtained by the Cambridge/ Aachen algorithm. The 
distribution from p-\-p collisions are scaled by Neinary [29 . Tri- 
angle symbols are from MB-Trig and corrected for efficiency, 
acceptance and energy resolution. Only statistical error bars 
are shown for the Au + Au data. Solid black squares are the 
distribution from p -\- p collisions, scaled by N Binary The sys- 
tematic uncertainty of the p -\- p jet spectra normalization is 
- 50%. 



4 Summary 

The full reconstruction of jets in 0-10% most central heavy 
ion collisions at RHIC energies is presented. Systemat- 
ics of the underlying heavy ion background subtraction 
is studied utilizing various algorithms and consideration 
of the jet area. The Neinary scaling is observed for the 
least-biased cuts with the given ~ 50% systematic uncer- 
tainty of the p + p jet spectrum measurement. An unbi- 
ased jet reconstruction in heavy ion collisions with MB- 
Trig data appears to be feasible. However, spectrum cor- 
rections are currently based on model calculations using 
PYTHIA fragmentation. This aspect, together with back- 
ground subtraction techniques, spectrum variations due to 
cuts and reconstruction algorithms, must be investigated 
further in order to assess the systematic uncertainties of 
this measurement. 

The first heavy ion run at LHC is expected in late 2009. 
The heavy ion background is predicted to be larger at the 
LHC than at RHIC, but there will be copious production 
of very energetic jets, well above background The 
large kinematic reach at the LHC may provide sufficient 
lever-arm to map out the QCD evolution of jet quenching 
[31j. The comparison of full jet measurements in the differ- 
ent physical systems generated at RHIC and the LHC will 
provide unique and crucial insights into our understanding 
of jet quenching and the nature of hot QCD matter. 
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Fig. 10. Jet yield per event vs Et for 0-10% central Au + Au 
collisions obtained by the LOHSC algorithm. The distribution 
from p+p collisions are scaled by Neinary |29j. Triangle symbols 
are from MB-Trig and corrected for efficiency, acceptance and 
energy resolution. Only statistical error bars are shown for the 
Au + Au data. Solid black squares are the distribution from 
p-\-p collisions, scaled by N Binary The systematic uncertainty 
of the p + p jet spectra normalization is ^ 50%. 



using PYTHIA as the fragmentation model, and/or may 
be an indication of modified fragmentation due to quench- 



ing. 
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